Dissipative tunneling and orthogonality catastrophe in molecular transistors 
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Transport through molecular devices with weak tunnel coupling to the leads but with strong cou- 
pling to a single vibrational mode is considered in the case where the vibration is damped by coupling 
to the environment. In particular, we investigate what influence the electrostatic coupling of the 
charge on the molecule to the vibrational modes of the environment has on the I-V characteristics. 
We find that, for comparable characteristic length scales of the van-der-Waals and electrostatic 
interaction of the molecule with the environmental vibrational modes, the I-V characteristics are 
qualitatively changed from what one would expect from orthogonality catastrophe and develop a 
steplike discontinuity at the onset of conduction. For the case of very different length scales, we 
recover dissipation consistent with modeling the electrostatic forces as an external influence on 
the system comprised of molecule and substrate, which implies the appearance of orthogonality 
catastrophe, in accord with previous results. 



I. INTRODUCTION 

The possibility of creating devices on a molecular level 
has opened up the field of single-molecule electronics in 
recent years. The transport properties of such meso- 
scopic systems have been investigated in a number of 
experimentsiSi^iiiSiSiLSi Of particular interest has been 
the influence of strong electron-phonon coupling on elec- 
tron transport,^'^'^ which manifests itself as emission and 
absorption of vibrational quanta observable in the excita- 
tion spectra. One example is the series of experiments by 
Park et al^ where it was shown that the current through 
a single Ceo molecule was strongly coupled to a single 
vibrational mode. 

A large amount of theoretical works has dealt with 
the problem of tunneling through a single molecular elec- 
tronic level that is coupled to phonon modes. Since in 
many experimental realizations the tunnel coupling of 
the molecule to the leads is rather weak compared to the 
other energy scales in the problem, transport is domi- 
nated by Coulomb-blockade effects and occurs sequen- 
tially. An approach based on rate equations can thus be 
justified and has been used in a number of recent pa- 
pgj.g|9ji0jiiji2 pjiygically, it is an essential question how 
the excited vibrational levels are allowed to relax, either 
through coupling to the environment, for example the 
phonons or plasmons of the metal substrate, or by virtue 
of the tunneling electronsiiii^^ In the case where the 
relaxation of the vibrational mode is faster than the tun- 
neling rate one can assume an equilibrium phonon distri- 
bution. 

The coupling between the vibrational mode of the 
molecule and the environment depends strongly on which 
vibrational mode is considered. For intra-molecular vi- 
brations the lifetime can be very long.^® -'^^ However, we 
were able to show in a recent paper^^ that the vibra- 
tional mode associated with the center of mass motion of 
the molecule is coupled more strongly to the environment 
and is thus exposed to an effective damping mechanism. 
The electron-phonon effects seen in the experiments of 



Ref. are suggested to be due to such a center of mass 
motion. 

In the present paper, we investigate how screening of 
the charge on the molecule influences the damping mech- 
anism. A sketch of the physical setup that serves as 
starting point is shown in Fig. ^ If the molecule is not 
charged, it is held in place by an interaction of the van- 
der-Waals type. It can perform center-of-mass oscilla- 
tions about its equilibrium position, and the character- 
istic energy huj of such an oscillation was found to be 
^ 5meVi^ These oscillations also influence the substrate 
to which the molecule is attached, and we will call this in- 
teraction the van-der- Waals interaction in the following. 
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FIG. 1: Illustration of the system considered in this paper, 
(a) The molecule is attached to a substrate, e.g. by van-der- 
Waals interactions. The molecule can perform center-of-mass 
oscillations about its equilibrium position, which also exerts 
forces on the substrate. We label this interaction the van- 
der-Waals interaction in the following, and we focus on the 
linear regime where the interaction strength is independent of 
the position of the molecule, (b) When an electron hops onto 
the molecule, the force created by induced surface charges 
causes a shift of the equilibrium position of the oscillator and 
of the vibrational modes of the substrate. For a pointlike 
charge on the molecule, the shift of the equilibrium positions 
takes place in an effective dipole fleld (multipole field for non- 
pointlike charges). We label this interaction the electrostatic 
interaction. 
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On the other hand, if the molecule is charged, there will 
be additional electric forces on the molecule, which are 
either due to static electric charges on the surface, e.g. 
because of impurities, or due to surface charges in the 
substrate that are induced by the charge on the molecule. 
Hence, when the electron tunnels onto the molecule, it 
results not only in a force on the molecule itself but also 
in electrostatic forces on the substrate to which it is at- 
tached, and the equilibrium positions of the vibrational 
modes in the substrate are shifted accordingly. In the 
following, we will focus on the case of surface charges 
only, but static impurity charges can be treated analo- 
gously and have the very same qualitative effect. We 
will refer to this type of interaction as the electrostatic 
interaction.^'* 

Under the influence of the Coulomb interaction, the 
surface is attracted to the molecule (and vice versa), and 
substrate atoms will be shifted away from their equi- 
librium positions until van-der-Waals forces and electro- 
static forces are in equilibrium. This shift, however, oc- 
curs in an effective dipole field (or multipole field for 
charges that are not pointlike) since the field between 
molecule and substrate due to the additional charge on 
the molecule is equivalent to the case without surface 
charges but an image charge located inside the substrate, 
cmp. Fig. 121 Only in the limit of a very short ranged 
van-der-Waals interaction compared to the separation of 
charge and image charge could we approximate the elec- 
tric field at the surface as constant, and thus monopolar. 
Such a mismatch of length scales is rather unlikely for 
a situation as depicted in Figs. ^ and |2Ia). However, a 
different experimental geometry as in Fig. Efb) can in- 
troduce monopolar effectsii^ 

The shift of the equilibrium positions at for each 
phonon mode k can be calculated by resorting to the dy- 
namical matrix D(k, k') which appears in the harmonic 
expansion of the Hamiltonian (see also Ref . ITgl) : 

iJharm = ^ akD(k, k')ak, (1) 
kk' 

If we denote the electrostatic force on the substrate as 
G(r) and choose a phonon eigenmode basis such that 
D(k, k') is diagonal, we can see that force equilibrium 
demands 

D(k)ak^^G(r)e^''- (2) 

r 

where the right hand side is the Fourier component of 
the force G acting on the plane wave phonon state k. 
The long-wavelength behavior can give rise to a version 
of orthogonality catastrophe, i.e. zero overlap between 
the ground states corresponding to the shifted and un- 
shifted equilibrium positions. For small wavevectors k, 
we know that D(fc) ~ k^, see Ref.l20l Expanding G in a 
multipole expansion, and expanding the exponential for 
small k, we can see that ak ^ k~^ for a dipole force since 
the volume integral of any multipole field is zero by def- 
inition,^-'^ whereas ak ~ fc~^ for a monopole field since 




FIG. 2; (a) Outside the substrate, the electric field between 
the charged molecule and the substrate is equivalent to the 
case of an image charge located inside the substrate and no 
surface charges. The shift of the equilibrium positions of the 
vibrational modes of the substrate can thus be understood to 
occur in an effective dipole field for pointlike charges (multi- 
pole field in general). Since the shift of the equilibrium po- 
sitions then has an integrable long-range part (see text), the 
small-frequency part of the phonon spectrum is not shifted, 
and thus there will be no orthogonality catastrophe associ- 
ated with the tunneling process. The situation is different 
in (b), where the molecule interacts with image charges in 
both electrodes. If the two electrodes can be considered as 
mechanically independent, there will be long-range compo- 
nents of the displacement fields (the second electrode can be 
thought of as an external infiuence), and the orthogonality is 
restored. This case was considered in our previous analysis 
Ref. [13 

this corresponds to the constant part in the multipole 
expansion. This difference in long-wavelength behavior 
gives rise to nonzero (dipole) or zero (monopole) over- 
lap between the ground states in three dimensions, as 
one can see from the following simple argument. Let 77 
be the ground-state overlap between a set of harmonic 
quantum oscillators and its shifted counterparts, then 




where = mkWk is the oscillator quantum with £^ ^ 
k~^ for small k. Converting the sum to an integration, 
we see that the question of zero or nonzero overlap is 
decided by both dimensionality and the small-A: behavior 
of the equilibrium shifts ak. In particular, 

77 ~ exp ( - / d/c k'^al/A] , (4) 
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V/UJQ 

FIG. 3: Sketch of the difference between dipolar (dashed line) 
and monopolar case (solid line) in the I-V characteristics. We 
depict the situation for £o=2a>o and Vl ~ V/2 — —Vr in the 
notation of Sees. HTl 11111 and II VI The dipolar curve features a 
finite step of size A/ at the onset of conduction, whereas the 
monopolar curve rises according to a power law. 



where d is the dimensionaUty, and A is a normaUzation 
factor. Thus, in one dimension, we encounter orthogo- 
naUty catastrophe for both monopolar and dipolar shifts, 
whereas in two and three dimensions we see that dipolar 
equilibrium shifts feature nonzero overlap! 

It is this effect that we investigate in the present paper 
in the context of how a single vibrational mode of the 
molecule dissipates energy into a phononic environment. 
Again referring to Fig. El if the range of the van-der- 
Waals interaction D is comparable with the dipole length 
scale (set by the separation of charge and image charge), 
we expect nonzero overlap between the two ground states 
in accord with a dipolar equilibrium shift, and we should 
thus see a step in the I-V characteristics, cmp. Fig. |31 
If, on the other hand, the van-der-Waals range is very 
small compared to the dipole length scale, we expect an 
orthogonality catastrophe corresponding to the approxi- 
mately constant (monopolar) force. The latter can also 
occur for particular geometries, cmp. Fig.j^l^b). The cur- 
rent then behaves according to a power law at the onset 
of conduction. 

Another way of looking at the same problem is in 
terms of translational invariance. Suppose that the van- 
der-Waals range D and the dipole length scale are com- 
parable to each other. In this case, we could think of 
drawing a box around the whole system which contains 
both the van-der-Waals and the electrostatic interaction. 
Hence, all forces are internal, and we have a translation- 
ally invariant Hamiltonian. Small D, on the other hand, 
corresponds to the case that, from the point of view of 
the van-der-Waals interaction, the electrostatic interac- 
tion is extremely extended and spread out. Therefore, 
the electrostatic force acts approximately as an external 
force as far as the van-der-Waals interaction is concerned 
and translational invariance is broken. This case can be 
recovered in a physical setup where the molecule is lo- 
cated asymmetrically in between two substrates, cmp. 
Fig-EIb). Then, image charges in the substrate farther 
away can approximately have the effect of an external 



force on the substrate closer to the molecule. 

The paper is organized as follows. The model Hamilto- 
nian is defined in Sec.^ and in Sec. IIIII we recapitulate 
the expression for the current derived on the basis of rate 
equations in Ref. 12. The function that describes the tun- 
neling density of states is then solved in the case of equal 
spread of van-der-Waals and electrostatic interaction in 
the presence of the dissipative environment in Sec. IIVI 
Finally, a discussion and summary of the results can be 
found in Sec. IVI 



II. MODEL HAMILTONIAN 

We consider a model of one single spin-degenerate 
molecular level coupled to two leads. The single level is 
coupled to the vibrational mode of the molecule through 
the charge on the molecule, as are the vibrational modes 
of the substrate. The coupling between the oscillator and 
the environment is included as a linear coupling to a bath 
of harmonic oscillators. Similar to the model studied in 
our previous paper, the model Hamiltonian then reads 



H = H 



LB. 



ath 



bath 



H 



DB 



with 



H 



LB ^ka'^kry,a'^kij,a 
ka, a—L.R 

^bath = 



^'O , -'^ 2 2 



H 



Sbath 




(5) 

(6a) 
(6b) 

(6c) 
(6d) 

(6e) 
(6f) 
(6g) 



k(7, a—L.B 



where the c^^ ^ , Cj^^ ^ and the ^ g^j-e creation and 
annihilation operators for the leads and the molecule, 
respectively, Xq is the oscillator degree of freedom, {xj} 
describes the set of environmental degrees of freedom and 
rrij, ujj their respective masses and frequencies, ^0 is the 
onsite energy, and U is comprised of the Coulomb inter- 
action and the energy cost for having double occupancy 
on the molecule. The coupling constants for the electron- 
oscillator interaction and the oscillator-bath interaction 
are A and {Pj}, respectively. The lead electron energies 
are given by 



^ka — ^ko 



fJ-a 



(7) 
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where /x^, is the chemical potential of lead a. The tun- 
neling amplitudes t}^a in the tunnehng term Hj- could 
in principle also depend on the oscillator position, which 
should be a small effect for the experimental realizations 
in Refs. and ^^"^ 

We can explicitly write down how the /3j are connected 
to the displacement u(v) perpendicular to the substrate 
if we assume that the force acting between the molecular 
vibrational mode and the substrate phonons is spread 
over the substrate according to a normalized distribution 
/(r), i.e. / dr/(r) = 1. If we model the interaction as a 
harmonic spring potential, we can write 



Given Eq. H12|l . translational invariance of the total 
Hamiltonian with respect to shifting both xq and u(v) by 
some constant displacement vector is ensured automati- 
cally if g[y) is normalized properly, i.e. J drf/(r) — 1. 

We now want to relate the coupling constants in the 
boson-bath coupling to the finite damping of the vibra- 
tional mode, which can be accomplished by studying the 
classical equations of motion. After removing the bath 
degrees of freedom, we end up with the following equa- 
tion of motion in the frequency domain, neglecting the 
coupling terms Hj^q and -ff^bath which will be removed 
by a unitary transformation below: 



-Hint = 2 "^0^0 



dr/(r)u(r) 



(8) 



Expanding the square, we can combine the term 
\mQiij\ [/ dr/(r)u(r)] with the part of the Hamiltonian 
that describes the "free" environmental bath and diag- 
onalize the combination by introducing new coordinates 
Xj. These are related to u(r) by a unitary transformation 



as 



(9) 



with appropriately chosen coefficient functions Cj(r). 
Defining 



1 dr/(r)c,(r) 



(10) 



we immediately see that we recover the structure of our 
original Hamiltonian H. The advantage of the form of 
-ffint is that its translational invariance with respect to 
shifting both and w(r) by the same constant displace- 
ment is immediately apparent since the distribution func- 
tion /(r) is normalized to 1. 

The difference to Ref . is the newly included coupling 
term 7?£)bath between the charge on the molecule and the 
bath of vibrational modes, 



Dbath 



(11) 



It is again instructive to rewrite this term as a function 
of u(r) and a normalized distribution function g(y) that 
characterizes the form and range of the electrostatic in- 
teraction between molecule and substrate: 



Dbath 



= -A 



(12) 



so that the electrostatic force on the substrate G(r) men- 
tioned in Sec. U has the form G(r) — Ag(r) c?^c?^. 
Since u(r) can be written as a linear superposition of the 
eigenmodes Eq. we can express the coupling 

constants \j in terms of the distribution function g(x') as 



(13) 



Xj^^X drg{r)cj{r). 



[cj^ -uj^- S{uj)]xo{uj) = 0, 



where we have defined 



= — 



1 

mo 



E 



/3? 



1 



(14) 



(15) 



which is complex in general and gives rise to frictional 
damping and a frequency shift of the bare frequency loq. 
In Ref. ha, we explicitly calculated S{u!) for the case of 
a molecule attached to a semi-infinite substrate. We will 
use these results throughout this paper, again assuming 
that the van-der-Waals force J- exerted by the molecule, 
cmp. Fig.^a), is directed perpendicular to the substrate 
surface and cylindrically symmetric. Also, we assume 
that both the van-der-Waals and the electrostatic forces 
only couple to the substrate surface, i.e., g{r),f{r) ~ 
S{x), where the x axis is perpendicular to the substrate 
surface. The van-der-Waals force T on the molecule is 
thus understood to be of the form 



Xo 



2'Kr f {r)ux{r)dr 



(16) 



where Ux{r) is the perpendicular displacement on the sur- 
face of the substrate, and r is the radial direction in the 
surface plane. For explicit numerical evaluations later in 
this paper, we will use /(r) = e~^/^ /2ttD^ just as in our 
previous paper, with D on the order of the width Dq of 
the molecule 

We now eliminate the coupling terms Hjj^ and iJnbath 
of the Hamiltonian ((S} by a unitary transformation sim- 
ilar to the one used in the independent boson model,— 
at the cost of introducing displacement operators in the 
tunneling term. However, since we have a somewhat 
more complicated system here due to the couplin g to the 
bosonic bath, the unitary transformation in Ref. |22 has 
to be generalized. We define the transformation 

H^SHS\ 5 = e-'^"^ A^p,e + J2p/r (17) 



where — J2a^l-^cr- Using that 

^0 = ^0 ~ ^'^di = ~ ^j'''^d: (18) 



5 



both linear coupling terms, Hj^g and -ffubath, cancel if 
we set 



- ^...2 ^ .,■',2 (19) 



m„[^2 + 5(0)]' 
The Hamiltonian then transforms into 

H = Hi^j^ + Hjj + + i?bath + -H^Sbath + Ht, (20) 

where 



where we defined 



A(.;) ^ - 



(29) 



If the two distribution functions /(r) and ^(r) deviate 
only slightly from each other, Cj is small compared to Pj, 
and we see that A <g; 1. Also note that the qualitative 
frequency dependence of A(ix') is exactly the same as of 



Ht= E Vo(cL.ae'^'^.+4e-'^Cfe,,„) (21) 



ki7^ a—L^R 



and 



Hd^ ^o^dlda + Un^^n^i^, eo^S.Q--\l. (22) 



Here we have defined 



mo[cj2 + 5(0)] 



^^0 = ^0 



(23) 



where 



i>j 



(24) 



Thus, the Coulomb repulsion U is modified by both the 
phonon mediated interaction and the electrostatic inter- 
action between molecule and substrate. 

It is important to note that for the particular choice 
/(r) = g(r), i.e. when the distribution functions of elec- 
trostatic and van-der-Waals interaction are equal, we ob- 
tain from Eqs. and ifT^ : 



A/3, ^ 



l.^Q. (25) 



Furthermore, we can quantify the deviation of /(r) from 
g{v) by defining 

e(r) = g{r) - /(r), = -moul j dre(r)cj(r), (26) 

where J dre(r) — 0. We can thus rewrite Aj as 



III. CURRENT FORMULA 

As mentioned in the Introduction, the most impor- 
tant assumption here is that the tunneling rate is much 
smaller than all other time scales, which means that we 
can assume the vibrational degrees of freedom and the 
Fermi seas in the two electrodes to be in equilibrium at 
all times. For simplicity, we consider only two charge 
states and therefore let U = oo, which leaves us with 
only three states playing a role in the rate equations: 
empty, and occupied by either spin up or down. Treating 
the tunneling term Ht as a perturbation and accordingly 
employing Fermi's Golden Rule for the tunneling rates, 
the current / through the molecule is given by 



/ = 2e 



10-*- 01 Ol-"^ 10 



(r^i + r?i) + 2(rfo + rf„)- 



(30) 



Here, we defined F^- as the tunneling rate for tunneling 
from the state with occupation i to the state occupation 
j through lead a (see, for example, Ref . l2^ : 



(31a) 



r?o=r„ / ^F(c.)nJeo+^), 
roi=r„ / ^F(-a>)(l-nJeo+^)), (31b) 



with 



F{lu)^ dtF{t), F{t) ^ {e'^'-'^e-'^), (32) 



where the operator A is defined in Eq. (|17ll . Furthermore, 
we work in the wide-band limit with bare tunneling rates 
r^=27r^^ \tka\'^S{^k), and we intoduced the Fermi dis- 
tributions of the two leads na{e)={e^^'^~'^^''^ + 1)~^. 



IV. INFLUENCE OF THE COUPLING TO THE 
ENVIRONMENT 



moUjQ 



This immediately leads to 



(27) 



i-^[l~Am,i, = ^V}Z\ (28) 



In presence of coupling to the environment, the eval- 
uation of the function F{t) in Eq. (|32|l is in principle 
straightforward since the Hamiltonian is quadratic in the 
oscillator and bath degrees of freedom. We obtain 

^^(t)=exp(B(t)-B(0)), B{t) ^ {A{t)A{0))o, (33) 
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where the operator A is defined in Eqs. (|17|l and (|19|l . 
The expectation value (■ • ■ }o is to be evaluated with re- 
spect to H without the tunneling term. We follow the 
same procedure as in our previous paper— and use the 
fluctuation-dissipation theorem, 



B{uj) = -2Im[S^(cj)](l + UBito)), 



(34) 



to express B{t) in terms of the corresponding retarded 
Green's function 



B^{t) = -iemA{t),Am)o. 



(35) 



Here, nB{cj) — (e^'^ — 1) ^ is the usual Bose function. 
We obtain 



- 1 



uj^ — lOq — S{uj) 
e^mo A{0)S{uj) - Oj^A{uj) 
'l-A(O) uj^ - - S{uj) 



[1-A(0)]2 

PtUq 



[A{OfSiuj) - 2A(0)A(u;)] 



E 



1 



[1 - A(0)]2 ^ mjUj'^ {lo + iriY - uj'^ 



, (36) 



where we have defined S{ui) = S{uj) — 15(0), A(w) = 
A{llj) — A(0), and the experimentally observable renor- 
malized frequency ujq = ujq + 5(0). 

If we choose a very wide-ranged distribution g{r), as 
compared to /(r), then the function e(r) will be close to 
— /(r), see Eq. This leads to ej sa —f3j, and upon 

insertion into Eq. (|28|l we recover the ex pres sions for £ 
and ij that we found previously in Ref. [IJ since then 
A(a;) « S{uj). In effect, the electrostatic interaction is 
thus so widespread compared to the van-der-Waals inter- 
action that it can be considered as an external influence 
on the system. In this case, we are justified to neglect the 
term i?_Dbath in our Hamiltonian altogether since trans- 
lational invariance is broken by this external force. It is 
a matter of straightforward algebra to show that B-^{uj) 
indeed reduces to the same expression as in Ref. 

If, on the other hand, the spread of the two interactions 
over the surface is comparable to each other, i.e. g{r) « 
/(r), we have the situation where Ej <^ Pj, and A <g; 
1 in Eq. The first term in Eq. ^ is of zeroth 

order in ej, the second term is of first order, and all the 
remaining terms are of second order. Furthermore, it is 
important to see that the imaginary part of all terms 
except of the last one are convergent for small lu since 
the small-frequency dependence of ImA(ci;) is identical 
to that of linS{uj) ~ lo. However, the last term is of 
second order in e^, so that to leading order, lmB^{uj) is 
convergent for small to if the distribution functions /(r) 
and g{r) feature comparable ranges. 

How small the ej really have to be in comparison to the 
f3j (i.e., how similar in range /(r) and g{r) have to be) 



in order for such a first-order expansion to make sense 
can be deduced from the following argument: from our 
previous analysisi^ we know that the divergent terms in 
lmB^{uj) for small ui behave like C/uj, where C is a nu- 
merical constant of second order in the ej. The diver- 
gence shows its predominant influence in the large-time 
behavior of F{t): 



F{t) = exp 



duj 
2^ 



(e- 



l)BiLo: 



(37) 



For large times t, the integral in the exponent diverges 
logarithmically, where the physical cutoff for the inte- 
gration is given by the smallest frequency of the phonon 
spectrum Wmin <S 'i'o-^^ In order for a flrst-order expan- 
sion in the ej to hold, we therefore have to demand that 
the following be true: 



In 



\C\ < 1, i.e. 



LUq 



< e 



1/|C| 



(38) 



Even if ej <^ [3j is enforced only moderately, this condi- 
tion is not very restrictive. We should therefore expect 
a step in the conductance even if /(r) and g(r) arc more 
than just "slightly" different! 

In the following, we will examine the particular case 
/(r) = g{r), i.e. \j = XPj/nioUjQ, see Eq. (^511 . for which 
the small-frequency behavior changes quite dramatically. 
(This is the case most different from our previous analysis 
in Ref. J^) Then we have ij = 0, and thus 



j'^ — UJq — S{uj) 



which leads to 



B{u;) 



Ag uj'^ll + nB{u!)]lmS{u!) 



^0 [cj2 _ ^2 _ Ke5(w)] - [Im5(cj)]2 



1 , (39) 



(40) 



Here, we have defined the experimental parameter 



1 



1 X^luq 

o 4 ' 



(41) 



which is a measure of how much the equilibrium position 
of the molecular vibrational mode is shifted compared to 
the oscillator quantum yJfi/mQUjQ. Note the appearance 
of both the renormalized and bare frequency in Eq. I|41|) 
and thus Eq. ll^ . 

One can then find F{t) from Eq. H37() . For large times 
t, the integral in the exponent is convergent and yields 
a nonzero finite number since ImiS(aj) oc uj for small fre- 
quencies in the case of coupling to a semi-infinite sub- 
strate— In particular, we obtain 



hm F{t) = exp ^ / — 
V^o J-oo 27r 



[1 + rig (oj)] lmS{uj) 



(42) 
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FIG. 4: Size of the discontinuity |A/| at eV = ±2eo (cmp. 
Fig.|5Ji calculated for a Ceo molecule on gold, with symmetric 
tunneling contacts Vl = (in units of the maximum current 
through the molecule /max = 4/]v/3) as a function of g and 
the particular choices D/Do = 0.5, 1, 2, oo (bottom to top). 
The large D / Dq limit coincides with the case of no vibrational 
coupling to the environment for which we have |A/| oc e~^. 
The dependence on D / Do is rather weak but still visible. 



where we defined 

g{oj) = [uj^ -^l- Re5(cj)] ^ - [Im5(tj)]2 (43) 

for the convenience of shorter formulae. Thus, we expect 
a i5-function at w = 0, which in turn results in a finite 
steplike discontinuity in the 1-V curves at zero tempera- 
ture according to Eqs. (jSnil^ . We define 



!(t) - F{t) - hm F{t), 

t—^oo 

and thus obtain at zero temperature: 



(44) 



F{lj) =27r f hm F{t)] 5{lj) + f{uj) (45) 



with 



4g up' \a\S(uS) 



2n 



( hm Fit)) 
(w - ^3 lmS{uj - C) 



(46) 



This implies that /(lu) oc lu'^ for small ui. For small we 
can actually approximate /(w) by the first term, since 
the other terms are of higher order in g. 

The resulting discontinuity A/ in the I-V character- 
istics calculated from Eq. H3U|) and measured in units of 
In = eTLTB,/ {Tl + r^) is given by 



— = 2sign(F) ( hm F{t) . 
In \t-*oo / 2 + Ll/Lr 



(47) 



We explicitly calculate AJ for the case of a Ceo molecule 
attached to a gold substrate, which is shown in Fig. 0] as 
a function of the experimental parameter g. The curves 
are parametrized by the spread D/Dq of the vibrational 
force over the surface of the substrate, see Eqs. (|15l I42|l . 



However, the dependence of the jump AJ on D/Dq is 
only weak, and the size is comparable to the step size 
encountered in the case without coupling to the envi- 
ronment, for which we had the same formula as above 
with limt^oo F{t) = e~^}'^ The discontinuity appears 
at the voltage eV = ±2eo for symmetric bias voltage 
Vl,r = i-V/2, but is not present for eo = for reasons 
of symmetry. 

Unless suppressed for reasons of device geometry, such 
a discontinuity at onset of conduction should be visible 
experimentally for weak tunneling into and out of the 
molecule even when the effect of finite tunneling rates F 
is taken into account. Following the derivations by Meir 
and Wingreen)2^ the differential conductance line shape 
is given by A(lu)T, where A{lli) is the spectral function of 
the molecule. Taking the ground state phonon overlap i] 
from Eq. it can be shown that^^ 



AiLj)T 



r2|^|6 



t^2+r2|77|4/4' 



(48) 



Thus, the height of the differential conductance step 
scales as \r]\'^, whereas the characteristic width of the 
step scales as FjTyp. For small enough F compared to the 
width of the first Franck-Condon step in the I-V curve, 
we therefore see that we recover a sharp step-like discon- 
tinuity in the I-V curve. We can estimate F by virtue of 
the maximum current through the molecule Imax ~ 

eF, 

which yields TiV ~ O.l/ieV since Imax ~ lOOpA for the 
experimental setup by Park et al^ The Franck-Condon 
steps, on the other hand, feature a source-drain voltage 
width on the order of mV in the I-V curves, see Ref.^ so 
that the discontinuity in the I-V curve will remain sharp 
even for finite tunneling. 



V. SUMMARY AND DISCUSSION 

We found that inclusion of the effects due to screening 
of the charge on the molecule can change the qualitative 
form of dissipation with which the molecular vibrational 
mode is faced. Such screening effects can be mediated 
by static surface charges on the substrate, or by surface 
charges that are induced in the substrate by the charge 
on the molecule. If the distribution functions that char- 
acterize the form of the electrostatic interaction and the 
van-der-Waals interaction of the molecule with the sub- 
strate are not too much different from each other, we 
encounter a step in the I-V curves at the onset of conduc- 
tion. Such a discontinuity should be observable experi- 
mentally in the weak tunneling limit even with smearing 
due to finite tunneling effects. This is a result qualita- 
tively different from treating the electrostatic forces as 
an external influence on the molecule-substrate system. 
The step size is dependent on the experimental parame- 
ter g that reflects the shift of the equilibrium position of 
the molecular mode compared to the oscillator quantum 
toqWo and features an exponential dependence on g sim- 
ilar to what one would obtain for the case of no coupling 
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to the environment. The exponent is weakly dependent 
on the spread D/Dq of the van-der-Waals coupUng be- 
tween molecular mode and substrate modes, and we re- 
cover the decoupled limit for D/Dq — cx) as expected 
from our previous results in Ref. 12. 

However, we also find that the results of our previ- 
ous analysis are not changed qualitatively (but can be 
changed quantitatively) if the electrostatic interaction is 
very long ranged compared to the van-der-Waals interac- 
tion, and we recover a power law for the current at onset 
of conduction. In particular, the limit — > taken in 
Ref. ^2 is explicitly seen to coincide with treating the 
electrostatic interaction as an external influence on the 
system that breaks translational invariance. 

The data in Ref. do not exhibit a step in the I-V 
characteristics and thus feature agreement with the case 
of an external force as chosen in Ref. It is quite pos- 



sible, however, that in other experimental realizations, 
van-der-Waals interaction and electrostatic interaction 
are of comparable range, without the geometry of the 
system suppressing the effects that we found to arise due 
to the absence of external forces. According to our calcu- 
lations, one would then indeed observe a significant step 
in the I-V curve at the onset of conduction. 
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